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ABSTRACT 

We present maps and radial profiles of the gas temperature in the nearby galaxy clusters A2199 and A496, 
which have the most accurate ASCA spectral data for all hot clusters. X-ray images, temperature maps, and the 
presence of moderate cooling flows indicate that these clusters are relaxed and therefore can provide reliable X-ray 
mass measurements under the assumption of hydrostatic equilibrium and thermal pressure support. The cluster 
average temperatures corrected for the presence of cooling flows are 4.8 ±0.2 keV and 4.7 ±0.2 keV (90% errors), 
respectively, which are 10% higher than the wide-beam single temperature fits. Outside the central cooling flow 
regions and within r k, 0.7 Mpc covered by ASCA, the radial temperature profiles are similar to those of the 
majority of nearby relaxed clusters. They are accurately described by polytropic models with 7 = 1.17 ±0.07 for 
A2199 and 7 = 1.24;!;y*j'j for A496. We use these polytropic models to derive accurate total mass profiles. Within 
r = 0.5 Mpc, which corresponds to a radius of overdensity 1000, riooo, for these clusters (estimated from our 
mass profiles), the total mass values are 1.45 ±0.15 x lO''*/?"' Mq and 1.55 ±0.15 x lO''*/!"' Mq. These values 
are 10% lower than those obtained assuming constant temperature. On the other hand, the values inside a gas 
core radius (0.07 -0.13 Mpc) are a factor of ;> 1.5 higher than the isothermal values. The gas mass fraction 
increases significantly with radius (by a factor of 3 between the X-ray core radius and riooo) and at riooo reaches 
similar values of 0.057 ± 0.005 Z!"^/^ and 0.056 ± 0.006/?"^/^ for the two clusters, respectively. Our measured 
total mass profiles within riooo are remarkably well approximated by the Navarro, Frenk, and White "universal" 
profile. Since A2199 and A496 are typical relaxed clusters, the above findings should be relevant for most such 
systems. In particular, the similarity of the temperature profiles in nearby clusters appears to reflect the underlying 
"universal" dark matter profile. The upward revision of the mass values at small radii for the observed temperature 
profile compared to those derived assuming isothermality will resolve most of the discrepancy between the X-ray 
and strong lensing mass estimates. 

Subject headings: cooling flows — dark matter — galaxies: clusters: individual (A496, A2199) — intergalactic 
medium — X-rays: galaxies 



L INTRODUCTION 

Under a reasonable, but as yet not directly tested, set of as- 
sumptions that the hot intracluster gas is supported by its own 
thermal pressure and is in hydrostatic equilibrium in the cluster 
gravitational well, one can determine the total mass of a clus- 
ter, including its dominant dark matter component (Bahcall & 
Sarazin 1977; Mathews 1978). Because clusters are the largest 
collapsed objects in the Universe, their mass values are of great 
importance for cosmology. The cluster mass function and its 
evolution with redshift constrain the spectrum of the cosmo- 
logical density fluctuations and the density parameter 51o (e.g.. 
Press & Schechter 1974; White & Rees 1978; Bahcall & Cen 
1992; Viana & Liddle 1996). If the cluster matter inventory is 
representative of the Universe as a whole, as is expected, then 
by measuring the cluster total and baryonic mass and compar- 
ing it to the predictions of primordial nucleosynthesis, one can 
constrain (White et al. 1993). Comparison of independent 
cluster mass estimates, for example, by X-ray and gravitational 
lensing (e.g., Bartelmann & Narayan 1995) methods, provide 
unique insights into cluster structure and physics. A discrep- 
ancy between the different estimates may indicate significant 
turbulence or nonthermal pressure in the intracluster gas (e.g., 
Loeb & Mao 1994), or the effect of line of sight projections. 

For an X-ray measurement of the cluster mass, one needs 



accurate radial profiles of the gas density and temperature, as 
well as confidence that the cluster is in hydrostatic equilib- 
rium. The gas density profile for a symmetric cluster can read- 
ily be obtained with an imaging instrument, such as Einstein or 
ROSAT. Obtaining temperature distributions has proven to be 
more problematic, especially for hotter, more massive clusters. 
ASCA (Tanaka, Inoue, & Holt 1994) now provides spatially re- 
solved temperature data for nearby hot clusters (e.g., Ikebe et 
al. 1997; Loewenstein 1997; Donnefly et al. 1998; Markevitch 
et al. 1998 [hereafter MFSV] and references therein), although 
their accuracy is still limited. Outside the central cooling flow 
regions, the temperature decreases with radius in most stud- 
ied clusters. For a few clusters with more accurate temperature 
profiles, accurate mass profiles were already obtained (e.g., for 
A2256 by Markevitch & Vikhlinin 1997, hereafter MV). 

MFSV found that gas temperature profiles of nearby sym- 
metric clusters outside the cooling flow regions are similar 
when scaled by the virial radius and average temperature. The 
gas density profiles also are rather similar (e.g., Jones & For- 
man 1984; Vikhlinin, Forman, & Jones 1999). This sug- 
gests that the underlying dark matter profiles are similar In- 
deed, analytical work and cosmological cluster simulations 
(e.g., Bertschinger 1985; Cole & Lacey 1996; Navarro, Frenk, 
& White 1995, 1997, hereafter NEW) predict that the dark mat- 
ter radial profiles of most clusters in equilibrium should be sim- 
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TABLE 1 
Cluster Parameters'' 
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A2199 
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2.24 X 10" 


4.4 ±0.2 


4.8 ±0.2 


0.65±0.11 


2.9 ±0.3 


0.161 ±0.014 


3.6±0.5 


10 


A496 


. . 249 


0.700 


1.02 X 10" 


4.3±0.2 


4.7±0.2 


0.47 ±0.10 


3.1±0.3 


0.158±0.017 


3.9 ±0.6 
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"All values are for h = 0.5. ''Best-fit values for the ROSAT PSPC brightness profile excluding the central r = 3'. The gas density value is an 
extrapolation of this /3-model to the center. '^Wide-beam single temperature fit. ''Emission-weighted temperature excluding cooling flows. 
^Involves extrapolation to the area not covered by the temperature profile. 



ilar in units of the virial radius. It is interesting to see whether 
their predicted "universal" dark matter profile agrees with the 
observations. 

For all but a few clusters in the MFSV sample, the temper- 
ature data have insufficient accuracy for such a test. We there- 
fore selected two additional typical, relaxed, but less distant 
clusters, A2199 (z = 0.030) and A496 (z = 0.033), for a more 
accurate temperature profile and mass derivation using ASCA. 
These clusters are very ordinary in their X-ray luminosities and 
temperatures (T ~ 4.5 keV) and, similarly to most clusters, 
have moderate cooling flows (170 and 95 Mq yr"', respectively; 
Peres et al. 1998). The presence of cooling flows is suggestive 
of a relaxed cluster, while at the same time these flows are not 
so strong as to prevent accurate resolved temperature measure- 
ments with ASCA (see MFSV). A subset of the data presented 
here (observations of the central regions) was already analyzed 
by Mushotzky et al. (1995). We have since obtained offset ob- 
servations, and include the ASCA PSF correction in our analy- 
sis. Below we use ASCA and ROSAT data on these two clusters 
to derive their total mass profiles. We use Ho=50 km s"' Mpc"' 
{h = 0.5); the error intervals are 90%. 

2. ROSAT" PSPC DATA 

To derive the gas density distribution, we use ROSAT PSPC 
data. The archival observations of A2199 and A496 were an- 
alyzed as prescribed by Snowden et al. (1994) and using S. 
Snowden's code. To optimize the signal to noise ratio, we 
used Snowden bands 5-7 that correspond to 0.7-2.0 keV. For 
A2199, two observations of the same field were combined. The 
radial brightness profiles were then fit with a /3-model Sxir) oc 
(1 + / al)~^^^2 plus a uniform X-ray background within a ra- 
dial range of 3' - 50'. The inner radius of 3' approximately 
corresponds to the cooUng radius for both clusters (e.g., Peres 
et al. 1998) and encompasses all of the X-ray brightness excess 
due to the moderate cooling flows in the cluster centers. The re- 
sulting parameters of the gas density profile are given in Table 
1 and are typical (e.g., Jones & Forman 1984). The /3-modeI 
values for A2199 are similar to the results of Siddiqui, Stewart, 
& Johnstone (1998) using the same data. 

3. A5CA DATA 

A2199 and A496 were each observed by ASCA with one cen- 
tral pointing and two different 14-15' offsets from the cluster 
centers. Such a configuration has been chosen to cover the clus- 
ter to a radius where the mean overdensity is 500 (rsoo), while 
at the same time keeping the cluster brightness peak within the 
ASCA field of view to avoid stray light contamination. Observ- 
ing the clusters at different positions in the focal plane also re- 
duces the ASCA systematic uncertainties that dominate in the 

''http://heasarc.gsfc.nasa.gov/docs/asca/abc/abc.html 



temperature estimates. The offset positions were chosen to 
avoid bright foreground sources and also to cover representa- 
tive regions of these sfightly elliptical clusters. 

After the standard data screening (ABC Guide''), useful GIS 
exposures for the A2199 central and offset pointings were 31 
ks, 19 ks, and 22 ks, and for A496, they were 37 ks, 24 ks, 
and 24 ks, respectively (the corresponding SIS exposures were 
about a factor of 0.8 of the GIS exposures). For the temperature 
fits, all pointings for both GIS and SIS were used simultane- 
ously; different pointings and instruments fitted separately give 
consistent results. To derive the spatial temperature distribu- 
tions, we used the method described in detail in MFSV and ref- 
erences therein. This method accounts for the ASCA PSF and 
assumes that outside the cooling flow regions, the ROSAT PSPC 
image provides an accurate description of the relative spatial 
distribution of the projected gas emission measure, after a cor- 
rection of the PSPC brightness for any gas temperature varia- 
tions. It should be mentioned here that a recent discovery of the 
possibly nonthermal EUV and soft X-ray (E < 0.2 keV) emis- 
sion in A2199 should not affect the latter assumption in any 
significant way, since we use a relatively hard (0.7-2.0 keV) 
PSPC band where this excess is absent (Lieu, Bonamente, & 
Mittaz 1999). The absorption column was assumed uniform at 
the Galactic values (Nh = 0.9 x 10^° cm"^ and 4.6 x 10^° cm"^ 
for A2199 and A496); for our £ > 1 .5 keV spectral fitting band, 
any expected variations are unimportant. 

The analysis method propagates aU known calibration and 
other systematic uncertainties, including those of the ASCA 
PSF, effective area, ROSAT and ASCA backgrounds etc., to the 
final temperature values. All reported confidence intervals are 
one-parameter 90% and are estimated by Monte-Carlo simula- 
tions. 

4. RESULTS 

4.1. Temperature Maps 

The resulting two-dimensional projected temperature maps 
are shown in Fig. 1, overlaid on the ROSAT images. We show 
only sectors in which the temperature is accurately constrained. 
The maps show no significant azimuthally asymmetric varia- 
tions, and together with the brightness contours suggest that 
these clusters are well relaxed. These maps may be contrasted 
to the similarly derived, but highly irregular, temperature maps 
of merging clusters, e.g., A754 (Henriksen & Markevitch 1996) 
and Cygnus-A and A3667 (Markevitch, Sarazin, & VikhHnin 
1999). In the central regions, the maps clearly show low tem- 
perature regions that correspond to the previously known cool- 
ing flows (e.g., Stewart et al. 1984; Edge, Stewart, & Fabian 
1992). 
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Fig. 1. — ASCA projected temperature maps (color) overlaid on the ROSAT PSPC brightness contours spaced by a factor of 2. Sectors in 
which the temperature was derived are numbered in upper panels and the temperatures with 90% errors are given in lower panels. Different 
colors correspond to significantly different temperatures. Vertical dotted lines separate regions belonging to different annuli. Dotted horizontal 
lines show average temperature inside the annulus. For the central cooling flow regions, a single-temperature fit is shown. A white circle in the 
map shows a point source excluded from the fit. 



4.2. Radial Temperature Profiles 

Figure 2 shows the cluster projected temperature profiles 
in five annuli. For the central radial bin, we used a model 
consisting of a thermal component and a cooling flow with 
the upper temperature tied to that of the thermal component, 
both with free normalizations. The figure also shows wide- 
beam, single-temperature fits (Te = 4.4 ± 0.2 keV and 4.3 ± 0.2 
keV for A2199 and A496, respectively) and emission-weighted 
average temperatures excluding the cooling flow component 
(Tx = 4.8 ± 0.2 keV and 4.7 ±0.2 keV). The latter are calculated 
from these temperature profiles as described in MFSV. The data 
indicate a higher temperature in the central cluster regions (out- 
side the cooling flows) compared to the average temperature, 
and a temperature decline with radius. This is similar to other 
clusters; in fact, when the profiles for A2199 and A496 are plot- 
ted in units of Tx and virial radius, they lie within the composite 
profile obtained by MFSV for other nearby, relatively symmet- 
ric clusters (Fig. 3). Such typical temperature profiles, together 
with the typical gas density profiles and the presence of cool- 
ing flows, make A2199 and A496 representative examples of 
relaxed clusters. 

Outside the central cooling flow bin, the profiles in Fig. 
2 are described remarkably well by a polytrope, Tgas oc pj^^ 
(both temperature profiles appear slightly more concave than 
the polytropic fits, which is probably nothing more than a coin- 
cidence; note that they differ in a similar way from the compos- 



ite profile in Fig. 3). Assuming the ROSAT-devived /3-models 
for pgas, we find 7 = 1.17 ±0.07 and 7 = 1.24!gj^ for A2199 
and A496, respectively. Regardless of whether this fact has any 
physical meaning or is purely fortuitous, it simplifies the to- 
tal mass derivation by providing a convenient functional form 
for the observed temperature profile. We will use it in the next 
section, but first note that for the mass derivation, one needs 
a real (three-dimensional) gas temperature profile as opposed 
to projected on the plane of the sky that we have obtained. 
We show in Appendix that as long as the gas density follows 
a /3-model and the temperature is proportional to a power of 
density, a projected polytropic temperature profile differs from 
the three-dimensional profile only by a normalization. For the 
best-fit (3 and 7 values for A2199 and A496, the projected tem- 
perature profiles are factors of 0.94 and 0.92 lower than the 
three-dimensional profiles, respectively. 

5. TOTAL MASS PROFILES 

For the total mass determination, we will take advantage 
of the fact that the temperature profiles can be described by 
a polytropic functional form. >From the hydrostatic equi- 
librium equation for a sphericaUy symmetric gas distribution 

Pgdfi(r) oc (1 +r^/fl^)"i^ and a temperature profile T oc Pgas', the 
total mass within a radius r = xa^ is given by 



M(O = 3.70x lO^^Mg 



0.60 T{r) a, 3f3-fx^ 



1 keV 1 Mpc 1 +x2 



(1) 
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Fig. 2. — Radial projected temperature profiles. Crosses are centered on the emission-weighted radii. Vertical errors are for 90% confi- 
dence; horizontal error bars show the boundaries of the annulus. Gray bands denote a continuous range of temperatures in a cooling flow, 
and the central cross corresponds to the upper (ambient) temperature of the cooling flow. Dotted and dashed horizontal lines show wide beam 
single-temperature fits and average emission-weighted temperatures excluding the cooling flow. Smooth lines show polytropic fits to the values 
outside the central bin, 7 = 1.17 and 1.24 for A2199 and A496, respectively. Values of riooo and rsoo are calculated from the mass profiles 
obtained from these temperature data (see Figs. 5 and 6 below). 



(see, e.g., Sarazin 1988). A polytropic temperature decline thus 
corresponds to the following correction to an isothermal mass 
estimate Mi.n' 
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Fig. 3. — Temperature profiles for A2199 and A496 (symbols) 
overlaid on the gray band representing a composite profile for a sam- 
ple of 19 nearby relatively symmetric clusters from MFS V (their Figs. 
7 and 8). For this comparison, the profiles are normalized by their 
cooling flow-corrected average temperatures Tx and plotted in units 
of ri8o estimated from Tx using the relation of Evrard et al. (1996). 
Cooling flow bins are not shown. The darker band corresponds to 
a scatter of best-fit temperature values of the MFSV sample and the 
lighter band covers most of their 90% intervals. 



M(r) _ T(r) 



(2) 



where T is the average temperature (see also Ettori & Fabian 
1999). To calculate the 90% confidence bands on mass pro- 
files (as well as the confidence intervals on the values of 7 
above), we have fitted the polytropic model to the same sim- 
ulated temperature values in those annuli that were used to 
calculate the temperature error bars (see MFSV). These fit- 
ted polytropic models were substituted into equations (1) and 
(2) above and 90% confidence intervals of the resulting val- 
ues were calculated at each radius. The resulting correction 
factor to the isothermal mass estimate is shown in Fig. 4 as 
a function of radius. The corresponding profiles of the total 
mass, M, and the gas mass fraction, /gas = Mgas/M, are shown 
in Fig. 5. The coiTesponding ratio of the mean total density 
within a given radius to the critical density at the cluster's red- 
shift [pc = 3H^{1 +z)V87rG] is shown as a function of radius 
in Fig. 6. For both clusters, our mass profiles correspond to 
'"1000 ~ 1.0 Mpc and rjoo ~ 1.3-1.4 Mpc (the latter involves 
extrapolation to a region not covered by the temperature pro- 
file, see Fig. 2). Masses and gas fractions at several interesting 
radii are also given in Table 1 . 

At r = 1 Mpc, Mushotzky et al. (1995) obtained mass es- 
timates of 2.55 X 10"^ Mq for A2199 and 3.05 x 10'"* Mq 
for A496. These estimates are close to ours, even though 
Mushotzky et al. did not apply the ASCA PSF correction in their 
analysis. From the galaxy velocity data, Girardi et al. (1998) 
obtained, at r = 1 Mpc, masses of 5.4+|;9 x 10'"* Mq for A2199 
and 3.5 ± 1.8 x 10'^ Mq for A496 (their 68% errors are mul- 
tiplied by 1.65 to obtain 90% intervals). These are consistent 
(although for A2199, only marginally) with our more accurate 
values. 
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Fig. 4. — Polytropic correction to the mass profiles derived under the isothermality assumption (see §g). Gray bands correspond to 90% 
temperature errors of the observed profiles shown in Fig. 2 (excluding the central cooling flow bins). The constant temperature is taken to be 
equal to the cooling flow-corrected values. X-ray core radii a, are from the ROSAT PSPC data (§0), and nooo and rsoo are same as in Fig. 2. 



6. DISCUSSION 
6. 1 . Uncertainty of Our Mass Values 

Although the polytropic model is a good representation of 
the observed temperature profiles, it does not necessarily cover 
all possibilities that could be consistent with the data. Since 
our best-fit model adequately represents the temperature and its 
gradient, the best-fit mass values should be unbiased (as long 
as the hydrostatic equilibrium assumption is valid). However, 
at the extremes of the confidence intervals, some other accept- 
able temperature models may result in slightly different mass 
profiles. Therefore, our confidence intervals on mass may be 
underestimated. A more exhaustive method of mass modeling 
would be to assume a certain functional form of the dark mat- 
ter profile and find parameter ranges consistent with the data 
(e.g., Hughes 1989; Henry, Briel, & Nulsen 1993; Loewenstein 
1994; MV; Nevalainen et al. 1999a,b). However, for the rela- 
tively high-quality data on A2199 and A496, we have chosen 
a simpler approach with the polytropic models, without giving 
undue importance to the formal error estimates. Indeed, given 
the present accuracy of the temperature and density profiles, 
the formal statistical uncertainties are already too small to be 
physically meaningful (e.g., MV). Hydrodynamic simulations 
suggest that systematic uncertainties in the method itself, such 
as the possible deviations from spherical symmetry and hydro- 
static equilibrium (in the form of significant gas bulk motions), 
can give rise to rms mass errors of about 15-30% (e.g., Evrard, 
Metzler, & Navarro 1996; Roettiger et al. 1996; see a more 
detailed discussion in §4.2 of MV). Those simulations included 
merging clusters in the statistical sample, so the relaxed clusters 
A2199 and A496 should have mass eiTors on the lower side of 
these estimates. 

Another source of systematic mass uncertainty is the pos- 

^NFW included only dark matter in their simulations and it is unclear whether 
subtract the gas mass from our total mass profiles (assuming, for example, the cun 
by the same functional form with slightly different parameter values. 



sible deviation of the measured electron temperature from the 
local mean plasma temperature that enters the hydrostatic equi- 
librium equation. Markevitch et al. (1996) proposed such 
nonequality as a possible explanation for an unusually sharp ob- 
served temperature gradient in A2163. Later theoretical work 
(Fox & Loeb 1997; Ettori & Fabian 1998; Chieze, AHmi, & 
Teyssier 1998; Takizawa 1998) concluded that for relaxed clus- 
ters, this effect should not be significant within the radial dis- 
tances presently accessible for accurate X-ray temperature mea- 
surements (about half the virial radius). Therefore, it is safe to 
assume that the mass values within riooo obtained in this paper 
are unaffected by this complication. 

Finally, at these low redshifts, the unknown cluster peculiar 
velocity may introduce a noticeable distance and mass error 
(e.g., a 1000 km s"' velocity would correspond to a 10% er- 
ror in the calculated mass). To summarize all of the above, the 
true uncertainty of the masses of A2199 and A496 is probably 
greater than our formal ±10% estimates and perhaps closer to 
20-25% (90% confidence at riooo), and is dominated by sys- 
tematics. 

6.2. The "Universal" Mass Profile 

NEW have found that radial density profiles of equilibrium 
clusters in their cosmological simulations can be approximated 
by a functional form p{r) oc (r / r^T^l + r / rs)~^ . This form is a 
very good description of our observed total mass profiles in the 
range of radii covered by the temperature data, as shown in Fig. 
5. Normalizations and scale radii r, of the NFW profiles were 
selected to fit the observed mass profiles. For A2199, r, = 0.18 
Mpc, and for A496, r^ = 0.36 Mpc.^ Extrapolating the best-fit 
NFW profiles to greater radii, we obtain the NEW's concentra- 
tion parameter c = raoo/r? of about 10 and 6 for the two clus- 

inclusion of gas would significantly change the shape of the mass profiles. If we 
ly favored value of h = 0.65), the resulting dark matter profiles are well described 
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Fig. 5. — Total enclosed mass profiles and the corresponding gas fraction profiles. All values are for h = 0.5 (the total mass scales as h ' 
and the gas fraction as h^^^^). The 90% confidence bands (gray) correspond to those in Fig. 4. Best-fit profiles of the Navarro, Frenk, & White 
form are shown as solid lines in the mass panels. In all panels, dashed lines show the profiles obtained assuming a constant temperature (at its 
average, cooling flow-corrected value). The riooo values correspond to the best-fit NFW profiles shown here (see also Fig. 6). 



ters, respectively. According to the NFW's simulations, c and 
the total mass within r2oo are strongly correlated for a given cos- 
mological model. Our c and M200 values agree well with those 
for several cosmological models considered by NFW, includ- 
ing CDMA (for that model, our observed masses correspond 
to M20()/M* w 5-6 as defined in NFW). An isothermal pro- 
file would imply a less concentrated dark matter distribution 
than that suggested by the NFW simulations (see also Makino, 
Sasaki, & Suto 1998). 

The outer regions of other clusters for which relatively ac- 
curate mass profiles were derived from the ASCA temperature 
profiles (e.g., A2256, MV; A3571, Nevalainen et al. 1999b) are 
consistent with the NFW profiles as well, although the con- 
straints are poorer Thus, the similarity of the gas temperature 
profiles found by MFS V for nearby relaxed clusters (outside the 



cooling flow regions) appears to be due to the underlying "uni- 
versal" dark matter profile of the NFW form. Note that we do 
not consider cooling flow regions due to their unknown temper- 
ature structure. As noted in MV and Nevalainen et al. (1999a), 
in those few relaxed clusters without cooling flows where the 
gas density profile exhibits a flat core all the way to the cen- 
ter (e.g., A401), the dark matter cannot have an NFW central 
cusp because the gas halo would be convectively unstable (see 
also Suto, Sasaki, & Makino 1998). On the other hand, it is 
likely that in clusters that do have central dark matter cusps, the 
corresponding dip of the gravitational potential causes the gas 
density peak and acts as a focus for a cooling flow. 

6.3. Mass - Temperature Scaling 
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Our mass values within our riooo are a factor of 1.6- 1.8 be- 
low the scaling relation between the total mass and emission- 
weighted average gas temperature derived from the simulations 
by Evrard et al. (1996). The same is true for other clusters 
(e.g., A2256, MV; A2029, Sarazin, Wise, & Mai-kevitch 1998; 
A401, A3571, Nevalainen et al. 1999a,b). Note that isother- 
mal mass estimates are also lower than the Evrard et al. M-T 
relation predicts. Given the agreement of our observed total 
(or dark) mass profile with the "universal" profile from NEW 
as well as from Evrard et al., the main source of this discrep- 
ancy apparently lies in the gas density and temperature distri- 
butions. Indeed, as noted by MFSV, simulations predict a less 
steep temperature decline than observed, and steeper gas den- 
sity profiles than observed (e.g., Vikhlinin et al. 1999). Both 
these effects have the sign needed to cause the M-T discrep- 
ancy. Because most of the cluster X-ray emission originates in 
the central region, simulations that do not sufficiently resolve 
the cluster core may predict an incorrect (apparently too low) 
emission-weighted gas temperature. For example, the Evrard et 
al. (1996) simulations have a resolution of 0.2-0.3 Mpc, com- 
parable to a typical cluster core radius. On the other hand, the 
hydrostatic mass measurement can underestimate the true total 
mass if, for example, there is significant gas turbulence. Indeed, 
simulations suggest that there may be residual turbulence even 
in an apparently relaxed cluster (e.g., Evrard et al. 1996; Nor- 
man & Bryan 1998) resulting in a 10-15% underestimate of the 
mass within r ^ riooo- 

6.4. Implications for the X-my-Lensing Mass Discrepancy 

A2199 and A496 are representative examples of relaxed clus- 
ters and the difference in their X-ray mass estimates using the 
measured temperature profile from the isothermal estimates, 
shown in Fig. 4, generally applies to other such clusters (see 
MV and MFSV). The upward revision of the mass estimate in 
the inner part has one important implication, the convergence 
of the X-ray and gravitational lensing mass estimates in the 
cluster central regions. The strong lensing mass values (that 



usually correspond to r ^ 0.2 Mpc) often exceed by a factor of 
2-3 the X-ray estimates made under the assumptions of isother- 
mality and a typical /3-model density profile (e.g., Loeb & Mao 
1994; Miralda-Escude & Babul 1995; Tyson & Fischer 1995). 
In many cases, the lensing analysis is likely to overestimate the 
mass as a result of substructure or projection (e.g., Bartelmann 
1995). On the X-ray side, some clusters are undergoing mergers 
and the hydrostatic equilibrium method may give a wrong mass 
value. For those clusters which are relaxed, the low-resolution 
X-ray image analysis may underestimate the gas density gradi- 
ent at small radii typical of the lensing measurements and may 
be responsible for part of the disagreement (e.g., Markevitch 
1997; Allen 1998). If the cluster has a strong cooling flow, the 
overall temperature can be significantly underestimated if no al- 
lowance for the cool component is made (Allen 1998), although 
for most clusters this correction is within ^ 20% (MFSV). 
Still, in many cases these effects alone are not sufficient to ac- 
count for the mass discrepancy. It has been suggested, e.g., by 
Miralda-Escude & Babul (1995) that a gas temperature gradient 
could explain the discrepancy for the distant non-cooling-flow 
cluster A2218. A temperature decline with radius has indeed 
been observed in A2218 by Loewenstein (1997) and Cannon, 
Ponman, & Hobbs (1999) using ASCA, while MFSV find that 
such a declining profile is common among nearby clusters. 

The analysis in §|| has shown that within the core radius, 
the commonly observed temperature gradient implies a mass 
that is higher than the isothermal estimate by a factor of ^ 1 .5. 
The reference isothermal estimate uses a cooling flow-corrected 
temperature, so this effect is in addition to the cooling flow- 
related mass correction of Allen (1998). If a similar temper- 
ature gradient is common in more distant clusters, this effect, 
together with others mentioned above, effectively resolves the 
mass discrepancy. This seems to obviate the need for more 
exotic causes, such as a significant magnetic field pressure or 
strong turbulence within cluster cores (Loeb & Mao 1994). 

6.5. Gas Fraction 



8 



MARKEVITCH ET AL. 



The lower panels in Fig. 5 show the gas mass fraction a func- 
tion of radius for the two clusters. At rjooo^ we obtain sim- 
ilar values of /gas = 0.161 ±0.014 and 0.158 ±0.017 for the 
two clusters, respectively. These values are consistent with 
those for A2256, 0.14 ±0.01 at 

1000 5 obtained by MV using an 
ASCA temperature profile, and for A401 (O.lS^g jJ^) and A3571 
(0.16;!;ogj) from Nevalainen et al. (1999a,b). Our values are 
also similar to the median values for large samples of clusters 
analyzed using the isothermal assumption: /gas = 0.168 from 
Ettori & Fabian (1999, scaled to riooo), and /gas = 0.160 from 
Mohr, Mathiesen, & Evrard (1999, for clusters cooler than 5 
keV). The latter similarity is due to the fact that the effect of 
the radial temperature decline on mass is small at this radius; 
at greater radii, the isothermal analysis underestimates the gas 
fraction as Fig. 5 shows. 

The values of the cluster gas fraction from X-ray analysis are 
often used to place constraints on the cosmological density pa- 
rameter (ilo < 0.3), under the assumption that /gas in clusters is 
representative of the Universe as a whole (White et al. 1993 and 
many later works). However, /gas increases with radius even if 
one assumes a constant gas temperature (e.g., David, Jones, & 
Forman 1995; Ettori & Fabian 1999), and the true increase is 
steeper as seen in clusters with measured temperature profiles. 
Fig. 5 shows that /gas increases by a factor of 3 between the 
X-ray core radius and riooo and does not show any evidence 
of flattening at large radii and hence asymptotically reaching 
a universal value. Although at some radius, the cluster must 
merge continuously into the infalling matter with a cosmic mix 
of components, that presumably happens at the infall shock ra- 
dius of ^ 2-3 riooo, well beyond the region presently accessible 
to accurate measurements. Note that both the dark matter and 
the gas mass within a given radius, and thus /gas, are dominated 
by the contribution at large radii. Cosmological simulations 
suggest that at smaller radii, a deviation from the universal /gas 
value is not large (e.g., Frenk et al. 1996). However, at the 



present stage, the simulations do not accurately reproduce the 
observed gas density, temperature and /gas profiles (e.g., MFSV; 
Vikhlinin et al. 1999). We conclude that our /gas values for 
A2199 and A496 are consistent with the constraints on flo de- 
rived in earlier works (e.g., Ettori & Fabian 1999; Mohr et al. 
1999), but caution that such estimates at present involve a large 
extrapolation. Future observatories Chandra and XMM will be 
capable of studying the cluster outermost regions and possibly 
determining the asymptotic value of /gas. 

7. SUMMARY 

The ASCA gas temperature maps and radial profiles for 
A2199 and A496 indicate that these systems are representa- 
tive examples of relaxed, moderately massive clusters. Our 
high quality temperature data imply total mass profiles that are 
in good agreement with the NEW simulated "universal" pro- 
file over the range of radii covered by the data (0.1 Mpc < r < 
/"looo ~ 1 Mpc). Because the temperature profiles of these two 
clusters are similar to the average profile for a large sample 
of nearby clusters in MFSV, this agreement indicates that the 
NEW profile is indeed common in nearby clusters. The upward 
revision of the total mass at small radii, by a factor of ^ 1 .5 
compared to an isothermal analysis, may reconcile X-ray and 
strong lensing mass estimates in distant clusters. The observed 
mass profile also implies a gas mass fraction profile steeply ris- 
ing with radius. While our /gas values at riooo support earlier 
upper limits on JIq, the steep increase of /gas with radius, not 
anticipated by most cluster simulations, suggests that we may 
not yet have correctly determined the universal baryon fraction 
and caution is needed in such analysis. 

This work was supported by NASA contracts and grants 
NAS8-39073, NAG5-3057, NAG5-4516, NAG5-8390, and by 
the Smithsonian Institution. 



APPENDIX 

PROJECTION OF THE POLYTROPIC GAS TEMPERATURE PROFILE 

Assuming a spherically symmetric gas density distribution of the form p{r) cx (1 + r^/a^)"i'^ (and taking a=l for clarity) and the 
polytropic temperature profile 

r(r)ap^-'(r), (Al) 

one can calculate a temperature profile that is emission-weighted (by p^) along the line of sight I, as a function of the projected 
distance from the center x (such that =x^ + /^), as 

7proj(^) = ^\oc cx // cx -^VL^^ , ' cx ^ \ = (1 +x2)-«7-i) ^ T(x). (A2) 

/o pHr)dl /o pHr)dl (l+x^ + P)-^l^ dl (l+xV^+J 



That is, the resulting projected temperature profile has the same shape as the real (three-dimensional) profile in eq. (|A1[). The relative 
normalization of the projected profile at jc = (and, therefore, at all radii) can easily be derived using the above formulae and is found 
to be: 

Tproj _r[|/3(l+7)-i] r(3/3) 



T r[|/3(l+7)] r(3/3-i)' 



(A3) 



The normalization of the projected temperature profile is sUghtly smaller than that of the three-dimensional profile. For (3 > 0.5 and 
7 < 5/3, their difference is less than ~ 20%. 

Strictly speaking, different temperatures along the line of sight are not simply weighted with p^; to obtain an exact projected 
temperature, a single-temperature fit to a multi-temperature spectrum should be performed in the ASCA energy band. However, 
the above similarity of the projected and real profiles holds for any weighting that is proportional to p^T", which approximates a 
wide range of possibilities. The normalization (A3) changes only weakly if ay^Q. For example, taking a = 0.5 (weighting with a 
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bolometric emissivity) instead of a = changes the normalization for our clusters by only 1 %, and by less than 5% for any reasonable 
(3 and 7. We have therefore assumed a = in §4-. 2 for simplicity. 
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